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1 Introduction 



Over past few decades there has been a lot of work on global analysis of 
nonlinear wave equations in n > 1 spatial dimensions [23], 

u tt = Au + f(x,u,Vu), u(t,x)eR, (WEa) 

iu t = Au + f(x, \u\, \Vu\)u, u(t,x)eC, (WEb) 

with one main focus being the study of blow-up phenomenon for the case of 
power nonlinearities / = ±|w| p , ±| Vu\ p . Some key tools with a number of uses 
in this study are conservation laws and symmetries. 

Conservation laws such as energy provide basic conserved quantities used in 
obtaining estimates on \u\ or |Vu| for smooth solutions, and also in defining 
suitable norms for weak solutions. Of considerable interest are extra conser- 
vation laws such as conformal energies that can appear for special powers p 
depending on the dimension n. Symmetries, in contrast, lead to exact group- 
invariant solutions and play a role in defining invariant Sobolev norms. Scaling 
symmetries are of special relevance, as the critical nonlinearity power for blow- 
up is typically singled out by scaling-invariance of a positive energy norm. 
Moreover, scaling transformation arguments give a means of relating the be- 
havior of solutions in different regimes, for instance, solutions at short times 
with large initial data can be scaled to long times with small initial data when 
the nonlinearity power is subcritical. 

In this paper we present classifications of conservation laws and symmetries 
admitted by a variety of physically and analytically interesting semilinear 
radial wave equations. The equations will be organized according to their 
variational structure as follows. 

First, we consider the standard (hyperbolic) nonlinear wave equation 

u tt = Au ± u p (NLW) 

as well as the nonlinear Schrodinger equation 

m t = Au ± \u\ p u (NLS) 

and its radial derivative variant 

iu t = Au ± i\u\ p (u r + r^u) (dNLS) 

V p + 2 J 

where A = r~ m d r r m d r = d^ + mr~ 1 d r is the radial Laplacian in n = m+ 1 > 1 
(spatial) dimensions. These PDEs each arise as the stationary points 5C/5u = 
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of a Lagrangian functional C = L[u]r m drdt, given by 

£nlwM = \{-u 2 t + u 2 r ) =F —^—u p+1 , 
2 p+1 

L>kls[u] = iuu t + \u r \ 2 =f — — |w| p+2 , 

p + 2 

^dNLsM = lUUt + | M r | 2 =F — |-^M P (mt r - UU r ) . 

p ~\~ A 

As is well known, the NLS equation also possesses a Hamiltonian formulation 
u t = i8TL/8u with = J °° H[u]r m dr given by 

H NLS [u] = -\u r \ 2 ±^-\u\ p+2 

p + 2 

where multiplication by i defines a Hamiltonian operator with respect to the 
L 2 Hermitian inner product on the radial line < r < oo. Next, while the 
derivative version (dNLS) of the NLS equation is not Hamiltonian, it does 
have a Hamiltonian variant 

hit = Au+ m( < m ~ 2 K ~ 2 u ± i ({\u\ p u) r + y M^r- 1 ^ (dNLS-H) 

which reduces to the standard derivative NLS equation in the case m = 0. 
This radial generalization arises from the Hamiltonian formulation 

Ut = r - m / 2 d r (r- m / 2 5H/5u), H = / °° H[u]r m dr (1) 

with . 

#dNLsM = ^(uU r ~ UU r ) ± — — l^l^ 2 - 

Here r - m / 2 Q rr - m / 2 i s easily verified to be a Hamiltonian operator with respect 
to the radial L 2 Hermitian inner product. In particular, such an operator T> 
with no dependence on u or u or their derivatives is Hamiltonian [20] iff 
it is skew-adjoint in this inner product so consequently the Poisson bracket 
associated to it by {V, Q}v = Jo° (SV/5u)V(5Q/5u)dr will be skew- Hermitian 
and obey the Jacobi identity for all real functionals V = / °° P[u]r m dr, Q = 
Io° Q[u]r m dr. Note that, as a consequence of the skew-adjoint property, the 
Hamiltonian 7i = J °° H[u]r m dr will formally be a conserved quantity, = 
(to within boundary terms at spatial infinity), for all formal solutions u. 

Last, based on the factorization of the Laplacian A = (d r + mr~ l )d r) we 
propose two radial generalizations of the modified Korteweg-de Vries equation, 

u t = {Au ± u p+1 ) r (mKdV-1) 

and 

Til 

u t = (Au ± u p+l ) r + —(Au± u p+l ) (mKdV-2) 

r 
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both of which have neither a Lagrangian nor Hamiltonian formulation except 
in the case m = 0. We also introduce a Hamiltonian variant given by 

u t = (Am ± u p+1 ) r + 2~(A« ± u p+1 ) (mKdV-H) 

with 

#mKdvM = ± — r^ uP+2 

using the previous Hamiltonian operator (JT]), specialized in the obvious way 
to real functions u — u. These radial mKdV equations are examples of third 
order evolutionary wave equations 

u t = AV x u + f(x,u, Vu,Att), u(t,x)eR, (WEc) 

with power nonlinearities / = u p, V x u, u p+1 , where V z = V is the radial 

gradient. Well-posedness of such wave equations is an interesting problem 
which has received some recent attention. 

We emphasize the wave equations (dNLS), (dNLS-H), (mKdV-l,2,H) for m ^ 
are new, being radial generalizations of the familiar (n — 1 dimensional) 
derivative Schrodinger equation and modified Korteweg-de Vries equation. 

To begin we recall the definitions of symmetries and conservation laws from 
an analytical perspective (see also [20][5] ). A point symmetry of a radial wave 
equation (WE) is a group of transformations given by an infinitesimal gener- 
ator 

5t = r(t,r,u), 8r = £(t,r,u), 5u = r](t,r,u) 
on the variables t, r, u in the real-valued case, or 

5t = r(t,r,u,u), 5r = £(t,r,u,u), 5u = r](t,r,u,u), 5u = fj(t,r,u,u) 

on the variables t, r, u and u in the complex- valued case, such that the equa- 
tion (WE) is preserved. On solutions, such a transformation in both cases is 
infinitesimally equivalent to 

5t = 8r = 0, 5u = 7] — ru t — <fu r (2) 

called the characteristic form of the point symmetry. The expressions rj, r, £ 
are determined by the Frechet derivative of the wave equation (WE) applied 
to 5u holding for all formal solutions u. (More precisely, one works in a jet 
space setting, using the coordinate space defined by t, r, u and all derivatives 
of u modulo the equation (WE) and its differential consequences. In jet space, 
a point symmetry is the prolongation of the operator X = rd t + £d r + rjd u . 
When u is complex-valued, jet space is enlarged in the obvious way by u and 
its derivatives.) The set of all infinitesimal point symmetries admitted by a 
given wave equation (WE) has the structure of a Lie algebra (for the operators 
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X under commutation). For a given point symmetry, invariant solutions are 
characterized by the form Su = r](t, r, u) — r(t, r, u)u t — £(t, r, u)u r = when 
u(t,r) G R, or Su = r](t,r,u,u) — r(t,r,u,u)u t — ^{t,r,u,u)u r = together 
with its complex conjugate when u(t,r) G C. 

A conservation law of a radial wave equation (WE) is given by a space-time 
divergence D t {r m '^> t ) + D r {r m ^> r ) that is equal to a linear combination of the 
equation and its differential consequences, so that 

Dt^ + D r ^ r + mr- l ^ r = (3) 

holds for all formal solutions u. The radial integral of the conserved density ^ 
formally satisfies 

d r°° 

— / ^r m dr = -r m # r |S° = 
dt Jo 10 

which vanishes when the flux \I/ r at spatial infinity is zero or decays faster 
than r~ m . Hence C = / °° v|/V m (ir formally yields a conserved quantity for 
the equation (WE). Conversely, any such conserved quantity arises from a 
conservation law ([3]). Two conservation laws are equivalent if their conserved 
densities differ by a radial divergence r~ m D r (r m Q) on all formal solutions u, 
giving the same conserved quantity C up to boundary terms. The set of all 
conservation laws (up to equivalence) admitted by a given wave equation (WE) 
forms a vector space, on which there is a natural action [9J by the Lie group 
of all admitted point symmetries of the equation (WE). 

Each conservation law ([3]) of a radial wave equation (WE) has an equivalent 
characteristic form where Dj(r m \l/*) + D r {r m ^> r ) is just proportional to the 
equation (WE) multiplied by an expression Q that depends on the jet variables 
(up to some finite differential order). Such expressions Q for which the product 
of equation (WE) and Q yields a total space-time divergence (and hence a 
conservation law on solutions of the equation (WE)) are called multipliers. 
There is a specific relation between multipliers Q and conserved densities 
in the case of a hyperbolic wave equation (WEa), conserved densities 
\&'(t, r, u, u r , Ut) modulo radial divergences correspond to multipliers 

Q(t,r,u,u r ,u t ) = 5(^ t r m )/Su t , (4) 

while in the case of an evolutionary wave equation (WEb) and (WEc), con- 
served densities ^{t, r, u, u r ) or \l/'(t, r, u, u, u r , u r ) modulo radial divergences 
correspond to multipliers 

Q(t, r, u, u r , u rT ) = 5(^V m ) /Su, (5) 
Q(t, r, u, u, u r , u r , u rr , u rr ) = — i5(\E f V m ) / Su, (6) 

respectively. In all cases the multiplier Q is determined by |2,3J the adjoint 
of the Frechet derivative of equation (WE) applied to Q, augmented by addi- 
tional equations formed from the Frechet derivative of Q itself, holding for all 
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formal solutions u. (More precisely, one works in the same jet space as for the 
computation of symmetries.) Thus the determination of conservation laws via 
multipliers is a kind of adjoint problem [I] of the determination of symmetries. 

If a wave equation (WE) possesses a Lagrangian formulation, then its multi- 
pliers Q define variational symmetries St = Sr = 0, Su = Q in the real- valued 
case or St = Sr = 0, Su = Q in the complex-valued case, such that the 
determining equations on Q reduce to conditions equivalent to those given 
by Noether's theorem [20|2|3] for the Lagrangian to be formally invariant (to 
within boundary terms at spatial infinity). Variational symmetries correspond- 
ing to multipliers (TjJ, (j5J), © for conserved densities containing derivatives 
of u or u will be a point symmetry Su = Q in the real-valued case only if 
the derivatives of u in Q at most are first order and appear linearly, or in the 
complex-valued case Su = Q will be a point symmetry only if the same is true 
for derivatives of u in Q and there are no derivatives of u in Q. 

For each of the wave equations (NLW), (NLS), (dNLS) and (dNLS-H), (mKdV- 
1,2,H), we now classify all admitted conserved quantities containing up to first 
order spatial derivatives, along with all admitted point symmetries. In these 
classifications, the nonlinearity power will be restricted to p 7^ 0, 1 for (NLW), 
and p ^ for (NLS), (dNLS), (dNLS-H), (mKdV-l,2,H), so that all cases 
of linear wave equations are excluded. No restrictions will be placed on m 
(even allowing non- integer values). All computations have been carried out 
using the computer algebra programs LiePDE and ConLaw [231 26] utilizing 
an enhanced version of the program CRACK [27] for solving overdetermined 
systems of equations containing parameters. 

Throughout we use the notation F(u, q) = j j uq ' q ^ . Also, a "±" sign will 
refer to the sign of the nonlinear term in the wave equations, usually called 
the focusing /defocusing cases, respectively. 



2 Point symmetries 

All of the wave equations (NLW), (NLS), (dNLS) and (dNLS-H), (mKdV- 
1,2,H) obviously admit time translations r = 1, £ = rj — 0, and, for m — 0, 
space translations ^ — 1, r — r] — 0. In addition the Schrodinger equations 
(NLS), (dNLS), (dNLS-H) admit phase rotations r = £ = 0, 77 = iu. 

The following tables list, firstly, the scaling symmetries admitted by these 
equations, and secondly, any extra admitted point symmetries. (Note these 
classifications exclude all linear cases i.e. p = 0; plus p = 1 for (NLW).) 
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Table 1. Scaling symmetries 



T 




V 


Equation 


t 


r 


2 

u 

p-1 


(NLW) 


2t 


r 


2 

u 

V 


(NLS) 


2t 


r 


1 

u 

V 


(dNLS), (dNLS-H) 


3t 


r 


2 

u 

P 


(mKdV-l,2,H) 



Table 2. Extra point symmetries for (NLW) 



r 




V 


Remarks 


r 


t 





Lorentz boost, m = 


t 2 + r 2 


2tr 


4 

tu 

p-1 


4 

inversion, m = 

p-1 



Table 3. Extra point symmetries for (NLS) 



T 




V 


Remarks 





2t 


iru 


Galilean boost, m = 


t 2 


tr 


( -r 2 - -t]u 


4 

inversion, m = — 1 H — 






V4 p J 


V 



The inversion and boost symmetries of equations (NLW) and (NLS) are well 
known [21], so our classification in tables [2] and [3] mainly provides a complete- 
ness result that these two equations do not possess any additional symmetries 
for special nonlinearity powers. 

Surprisingly, no extra symmetries are found to be admitted by equations (dNLS) 
and (dNLS-H). 

Table 4. Extra point symmetries for (mKdV-1) 



r 




V 


Remarks 





2t 




Galilean boost, m = 0, p = 1 



Table 5. Extra point symmetries for (mKdV-2) 



T 




V 


Remarks 





2t 




Galilean boost, m = 0, p = 1 


t 2 


2 

-tr 

3 


--(4tu ± r) 


inversion, m = 1, p = 1 



Table 6. Extra point symmetries for (mKdV-H) 



T 






Remarks 







Tl 


Galilean boost, m = 0, p = 1 


t 2 


2 

-tr 

3 


-g(4tu ± r) 


inversion, m = 2, p = 1 
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Like the (NLS) equation, the mKdV equations (mKdV-2,H) possess boosts 
and inversions as extra symmetries for special nonlinearity powers, while no 
extra symmetries are found to be admitted by equation (mKdV-1). Our results 
in tables HJ El El for these radial mKdV equations with m > are new. 

Note the powers for which the inversions exist are called a conformal power. 



Table 7. Conformal powers 



V 


Equation 


4 

1+ — 

m 


(NLW), m ^ 


4 

rn + 1 


(NLS), 


1 


(mKdV-2), m = 1; (mKdV-H), m = 2 



As a summary of the main results, we list the point symmetry algebras found 
for m > 0, with their generators denoted by X trans , X scal , X inver , X phase . 
These algebras fall into four classes (with a suitable normalization of X sca i 
in each case): For equation (mKdV-1) as well as the non-conformal case of 
equations (NLW) and (mKdV-2,H), the admitted algebra is generated by the 
time-translation and scaling symmetries 

[Xt rans , -Xs Ca j] Xi rans . 

In the conformal case of equations (NLW) and (mKdV-2,H) this algebra is 
enlarged by an inversion symmetry 

trans ) A; nver ] 2X SCSL \, \X sca \, X[ nver j X[ nver . 

For the conformal (NLS) equation, its admitted algebra is a central extension 
of the algebra for the conformal (NLW) equation, as generated by the phase 
rotation symmetry 

phase; At rans ] [Xphase > -^scal] [Aphagg, ^mver] 0) 

while the admitted algebra for the non-conformal (NLS) equation as well as 
for the (dNLS) and (dNLS-H) equations is a similar central extension of the 
algebra for the non-conformal (NLW) equation. 
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Table 8. Point symmetry groups 







TV mi ti 1" i / \ n 


-X"trans j -^"scal 


(7(1) xi U(l) solvable 


(mKdV-1), non-conformal (NLW), 
non-conforrneil ^mKd\^-2 ,H) 


-^"trans; -^scal) -^invcr 


SX(2,R) semisimple 


conformal (NLW), 
conformal (mKdV-2,H) 


^trans i ^scab -^phase 


(17(1) x 17(1)) x 17(1) 
solvable 


non-conformal (NLS), 
(dNLS), (dNLS-H) 


-^"trans i -^~scal j ^inven -^phase 


SL(2,R) x (7(1) centrally 
extended semisimple 


conformal (NLS) 



3 Optimal symmetry groups for invariant solutions 

All one-dimensional point symmetry subgroups Qn\ determine corresponding 
group-invariant solutions u(t, r) when the above wave equations are augmented 
by the invariant surface condition Xu = where X is the symmetry generator 
of in characteristic form. In particular, the invariant surface condition re- 
duces a wave equation (WE) to an ODE [20. 5J. Any two conjugate subgroups 
will give rise to reduced ODEs that are related by a conjugacy transformation 
in the full point symmetry group Q acting on the invariant solutions u(t, r) de- 
termined by each subgroup. Hence, up to the action of Q, all invariant solutions 
for a given wave equation can be obtained by selecting a one-dimensional sub- 
group in each conjugacy class of all admitted one- dimensional point symmetry 
subgroups £?(!). Such a selection is called an optimal set of subgroups [2"Tj . 

Optimal subalgebras have been classified in the work presented in [22] for all 
three- and four- dimensional Lie algebras. Applying this classification to the 
point symmetry algebras of the wave equations under consideration, we have 
the following list of optimal one-dimensional subalgebras for finding all invari- 
ant solutions of equations (NLW), (NLS), (dNLS) and (dNLS-H), (mKdV- 
1,2,H) for m > 0: 



Table 9. Optimal one-dimensional subalgebras 
(a, b, c are arbitrary constants) 



Generators 


Equation 


-Xtrans , -^scal 


(mKdV-1), non-conformal (NLW), 
non-conformal (mKdV-2,H) 


-^trans 5 ^scal j -^invcr 4* -^trans 


conformal (NLW), conformal (mKdV-2,H) 


-^phase j ^"scal + tt ^"phasc , -^trans + bX sca \ + cYphasc 


non-conformal (NLS), (dNLS), (dNLS-H) 


-^phasc 3 -^trans , -<^trans i -^"scal 3 

-^scal + a ^phase, Atrans + -^invcr + ^^phasc 


conformal (NLS) 
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Invariant solutions for (NLW), (NLS), (dNLS) with m > (i.e. in all dimen- 
sions n > 1) have been derived in [4"|23II16II11|I12|I10II17] . We will present invari- 
ant solutions of (dNLS-H), (mKdV-l,2,H) for m > (i.e. in dimensions n > 1) 
elsewhere. (The m = case of (mKdV-l,2,H) is treated in [2"0~1I18|I14||2"3"] ). 



4 Local conservation laws 

As the wave equations (NLW), (NLS) and (dNLS) each have a Lagrangian for- 
mulation, all their admitted variational point symmetries yield corresponding 
conserved quantities which are well known [21] for (NLW) and (NLS). In par- 
ticular, time translation symmetry yields energy, space translation and boost 
symmetries yield momenta, and inversion symmetry yields conformal energy, 
while phase rotation symmetry yields charge. In addition, there is a special 
nonlinearity power for which the scaling symmetry becomes variational and 
yields a dilational energy. Our classification in the following tables provides a 
completeness result that no additional conserved quantities up to first order 
are admitted by these equations (NLW), (NLS) and (dNLS) for special non- 
linearity powers (excluding all linear cases i.e. p = 0; plus p = 1 for (NLW)). 

The remaining wave equations (dNLS-H) and (mKdV-l,2,H) do not have a La- 
grangian formulation, so consequently their admitted conservation laws come 
from multipliers given by adjoint-symmetries [3] rather than symmetries. For 
the Hamiltonian equations (dNLS-H) and (mKdV-H), note the Hamiltonian 
itself provides one conservation law. Our results in the following tables for all 
these radial equations with m > are new. 



Table 10. Local conservation laws for (NLW) 







Remarks 




— u t u r 


energy 


U t U r 


-i(«2+«2) T F(«,p + l) 


momentum, m = 


~r(u 2 + u 2 ,) + tu t u r + rF(u,p+ 1) 


-^t(u 2 + u 2 ) - ru t u r 
=FtF(u,p+l) 


boost momentum, 
m = 


-t(u 2 + u r ) + ru t u r 
2 2 

H uu t + tF(u, p+1) 

P- 1 


r(u 2 + u 2 .) — tutu r 

2 2 

uu r =p rF(u, p+1) 

p-1 


dilational energy, 
4 

m = 

p-1 


-(t 2 + r 2 )(u? + u 2 ) + 2tru t u r u 2 

2 p — 1 

H tuu t T U 2 + r 2 ) F{u,p + 1) 

p — 1 v ' 


-tr(u 2 + u 2 ) - (t 2 + r 2 )u t u r 
4 

tuu r =p 2trF{u, p+1) 

p-1 


conformal energy, 
4 

m = 

p-1 
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Table 11. Local conservation laws for (NLS) 





\5,r 


Remarks 




l(u r U — UU r ) 


charge 


iK| 2 T F(H, P + 2) 


~{u r Ut - U t U r ) 


energy 


it(u r u — u« r ) — r|it| 2 


\t(uut — UtU) + \r(uu r — U r U) 

i nil l2 1 a j. 7~i / 1 1 i r»\ 

+2t\u r \ ± itp (\u\,p + 2) 


boost momentum, m = 


^(u r u — uu r ) 


-(uut — utu) 
+K| 2 ±2F(|u|,p + 2) 


momentum, m = 


1 [ 

-t\u r \ 2 H — r(u r u — uu r ) 

2 8 

=F*F(M,p + 2) 


t(u t U r + U r Ut) H«r| 2 

2 4 

+ — r(utu — uut) (uu r + u r u) 

8 2p 

T -rF(\u\,p + 2) 


4 

dilational energy, m = 1 

P 


i t 2| Ur |2_ l r 2 M 2 

+ — tr(u r u — uu r ) 
T t 2 F(\u\,p + 2) 


— — t 2 (utU r + U r Ut) — —t\u r \ 2 

+ — r 2 (u r u — uu r ) H — triutu — uut) 
8 4 

t(u r S + im r ) =(= t7\F(|u|,p + 2) 

2p 


conformal energy, m = 1 

P 



Table 12. Local conservation laws for (dNLS) 







Remarks 


\u\ 2 


m r « — iuu r =F 2-F(|u|,p + 2) 


charge 


\u r \ 2 =F i(u r u - «n r )|«|~ 2 F(|n|,p + 2) 


— (u r Ut + Mt«r) 

=Fi(«t«- uu t )\u\- 2 F(\u\,p + 2) 


energy 


-(u r U — UU r ) 


-(UUt — Utu) + |«r| 2 


m = 


— {u r U — UU r ) + 2t\u r \ 2 

^2it(u r u - uu r )\u\~ 2 F(\u\,p + 2) 


— (UtU — UUt) — 2t{u r Ut + UtU r ) 

±2it(u t u - uu t )\u\-' 2 F(\u\,p + 2) 

(u r U + UU r ) — r\u r \ 2 

p 


dilational energy, 
2 

m = 1 

P 



Table 13. Local conservation laws for (dNLS-H) 







Remarks 


r m / 2 (u + u) 


r- m / 2 (i(« r - it r ) =F \u\ p (u + «) + — r" 1 ^ - «)) 
2 




r- m / 2 i(u- it) 


r -m/2 ( - nr + u r ± l\u\P{u -U) + (u + u)) 




— (uu r — U r U) 

±2F(|u|,p + 2) 


— (UtU — UUt) — Wr\ 2 ~ r ~ 2 \ u \ 2 

— =p i|w| p (wn7 — u r u) 


Hamiltonian 


\u\ 2 


\{uu r - u r u) =F 2(p + l)F(\u\,p + 2) 


charge, m = 


— tiuUr — U r U) 
2 

±2tF{\u\,p + 2)+ r -\u\ 2 


^t(u t u - uu t ) - t(\u r \ 2 + ^-r-' 2 \u\ 2 ) 

-t(\u\P+ 2 ± i|«|P(MU r - U r «)) 

+ ^r(u r u - «u r ) T (P + l)rF(|u|,p + 2) 



2 

m = 2 

P 
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Table 14. Local conservation laws for (mKdV-H), m ^ 







Remarks 


r -m/2 u 


_ r -m/2 ^ + mr -i Ur ± u p+i) 


mass 




1 9 1 

—u t u r + —u^. r + rnr L u r u rr 

TTl 2 L 
±U p+1 U r r H r" 2 «2 ± mr -l„P+l„ r + - u 2 P+2 

2 2 


Hamiltonian 


3 1 

-tul ru 2 =F 3F(u,p + 2) 

2 2 v y 


-3tu t U r + ^tu 2 r + ^tu 2p+2 ± 3ttt p + 1 Mrr + rUMrr 

p - 2 , (p - 2) 2 „ „ 1 , 

-6- tr- 1 u r u rr + 6 „ fr~ 2 < ru? 

p p 2 2 

-2^— UU r ± (p + l)rF(u, p + 2) 
P 


4 

m = 2 

V 




y 

-9t 2 r _1 u r Ui + -t 2 r _1 u 2 r + 18t 2 r~ 2 u,.«rr ± 3r _1 « 




-t 2 u 2 =F 3t 2 « 3 - 3tru 2 T 
2 

^f6tr~ 1 u =p r 2 u 


±9t 2 r _1 U 2 ltrr + 6tMMrr ± Ur^Urr + 6tr~ 2 U 2 

+18t 2 r~' i u 2 - 3tu 2 ± 18t 2 r- 2 u 2 u r + r 2 u 2 
g 

±l2tr~ 2 u r T u r + -t 2 r -1 M 4 ± 4tu 3 ± ru rr 
2l 


m = 2, p = 1 



Table 15. Local conservation laws for (mKdV-1), m/0 







Remarks 




—r~ m (u rr + m,r _1 tt r ± 


mass 



Table 16. Local conservation laws for (mKdV-2), m^O 







Remarks 


r~ m u 


— r~ m (u rr + rar _1 u r rfc u p+1 ) 


mass 


l r ^/ 2 u 2 
2 


-r 1/2 uu rr + ir 1 / 2 « 2 - r _1 / 2 uu r T -r 1 / 2 ^ 4 ± i r ~ 3 / 2 u 2 
2 4 8 


dilational momentum, 
m = 3/2, p = 2 


1 2 
-rtr 

2 


1 9 3 . 
—ruu rr + ~ ru r ~ 2uu r =p —ru 


dilational momentum, 
m = 3, p = 2 



For completeness we mention that in the m = case the well known mKdV 
conservation laws are listed in 119,18,20,21. 



5 Norms and critical powers 

The wave equations (NLW), (NLS), (dNLS) and (dNLS-H), (mKdV-H) possess 
both a scaling symmetry (which is uniform in m) and a conserved energy or 
Hamiltonian, so they each have an associated critical power with respect to 
the energy norm, E[u] = e[u]r m dr, as shown in the following table. 
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Table 17. Energy norms 



E[u] 


critical power p 


Remarks 


Jo ( 2 (|nr ' 2 + ' Ut|2) T F(U ' P + 1} J rmdr 


m — 1 


(NLW), m ^ 1 




4 

m — 1 


(NLS), m ^ 1 


XT (i |Ur|2 T " ™ u r)\u\- 2 F(\u\,p + 2))r m dr 


2 

m — 1 


(dNLS), m ^ 1 


J °°Q(un r - uu r ) ± 2F(\u\,p + 2)jr m dr 


2 

m 


(dNLS-H), m ^ 


J OO Q|« r | 2 ±F(n,p + 2))r"'dr 


4 

m — 1 


(mKdV-H), m/1 



All these wave equations also possess dilational energies or dilational Hamil- 
tonians for special nonlinearity powers p depending on m, in addition to the 
well known conformal energies admitted for the (NLW) and (NLS) equations 
in the case of conformal powers p (cf. table [7]). Interestingly, the Hamiltonian 
mKdV equation (mKdV-H) admits a conformal energy in this case too. 



Table 18. Dilational energies and dilation powers 



Dilational energy 


Dilation power p 


Remarks 


I te[u] + (ru r + 




r m dr 


4 

1 + — 

m 


(NLW), m ^ 


f °° ( te[u] H r(uu r — u r u) 

JO \ 4 J 


r m dr 


4 

m + 1 


(NLS), 


r°° ( te\u] H r(uu r — u r u) 

jo y 4 ( 


r m dr 


2 

m + 1 


(dNLS), m^-l 






2 

m + 2 


(dNLS-H), m ^ -2 


/ °° (*eH - ir« 2 ) 


r m dr 


4 

m + 2 


(mKdV-H), m ^ -2 



Table 19. Conformal energies and conformal powers 



Conformal energy 


Conformal power p 


Remarks 


Jo" 


(i 2 + r 2 )e[it] +2i(r« r + -^u)u t jr m dr 


4 

1 + — 

m 


(NLW), m ^ 


JT 


t 2 e[u] + ~tr(u-u r — «r«) — ^»* 2 M 2 ^ 


r m dr 


4 

m + 1 


(NLS), 


Jo" 


t 2 e[u] — i (tru 2 + itr _1 u + ^ r2u ) 


j r 2 dr 


1 


(mKdV-H), m = 2 



An interesting pattern in tables 18 and 19 is that when m is expressed in 
terms of p then the difference m cr it. — man. (for any fixed p) is equal to 
ord(<9 r ) — ord(<9() + 1 > where "ord" refers to the highest order of a speci- 
fied derivative appearing in the wave equation. Accordingly, the dilation and 
conformal powers are subcritical in all cases. 

Two other norms of analytical interest are the radial L 2 norm and the radial 
H s norm given by ||w||i,a = (Jq° \u\ 2 r m dr) 1 ^ 2 and = (Jq 00 |9*M| 2 r m <ir) 1 ^ 2 

for any positive integer s. The latter norm has a natural extension to all s > 
defined in terms of the Fourier transform u = / K m+i u(t, \x\) exp(—k ■ x)d m+l x 
such that u is in H s iff (1 + \k\ s )u is in L 2 (R m+1 ). The following tables list 
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the critical powers p for which these norms are scaling- invariant. Note in the 
case of the Schrodinger equations (NLS), (dNLS) and (dNLS-H), the L 2 norm 
coincides with the conserved charge. 



Table 20. L 2 critical powers 



critical power p 


Remarks 


1+ \ 
m + 1 


(NLW) 


4 

m + 1 


(NLS), (mKdV-l,2,H) 


2 

m + 1 


(dNLS), (dNLS-H) 



Table 21. H s critical powers 



critical power p 


critical s 


Remarks 


4 

1 + 


m + 1 2 


(NLW) 


m + 1 — 2s 


2 p-1 


4 

m + 1 - 2s 


m + 1 2 
2 p 


(NLS), (mKdV-l,2,H) 


2 

m + 1 - 2s 


m + 1 1 
2 p 


(dNLS), (dNLS-H) 



6 Concluding remarks 

The utility of symmetries and conservation laws can be extended by means of 
potential systems [EKE]. A potential system for a radial wave equation (WE) 
arises from any conservation law such that vanishing set of its multiplier, 
Q = 0, is contained in the set of all formal solutions u of the given equation. 
Potentiating such a conservation law yields the system 

vt = r m ^ r , v r = -r m ¥ 

whose solutions v up to shifts (v — > v+c for an arbitrary constant c) are in one- 
to-one correspondence with the set of solutions u. For a given potential system, 
any admitted symmetry or conservation law that has an essential dependence 
on the potential v represents a nonlocal symmetry or nonlocal conservation 
law, respectively, of the wave equation (WE). Of course, the superposition of 
a local symmetry or a local conservation law with a nonlocal one yields further 
nonlocal ones, and so for the purpose of classifications we will mod out the 
admitted sets of local symmetries and conservation laws. 

All potential systems arising from the conservation laws for the wave equations 
in tables [10] to [16] are given by potentiating: the mKdV equations themselves 
(mKdV-l,2,H) and the Hamiltonian variant of the derivative Schrodinger 
equation itself (dNLS-H); the charge conservation law for the Schrodinger 
equations (NLS), (dNLS), and (dNLS-H) in the case m = 0; and the dila- 
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tional momentum conservation laws for the mKdV equation (mKdV-2) in the 
cases m = 3, |. We find that only the first of these potential systems — poten- 
tiation of the (mKdV-1) equation itself (including the case m = 0) — yields 
nonlocal conservation laws and none yield any nonlocal symmetries. 



Table 22. Nonlocal conservation laws for (mKdV-1) 
e = -\/=fT = i, 1 respectively in the defocusing/focusing cases 



Potential system 


-Jri 




Remarks 


v r = u, 

Vt = U rr — S 2 U S 


e <ss/2v 


-e e ^(u 2 - V2u r ) 


m = 0, p = 2 


V r = U, 

Vt = Urr H r _1 « r — £ 2 « 3 

2. 


re sV2v 




m = 3/2, p = 2 


V r = U, 

Vt = Urr + 3r~ 1 u r — e 2 u 3 






m = 3/2, p = 2 



Each of these nonlocal conservation laws gives rise to a further potential sys- 
tem. Of the three, none yield any additional nonlocal conservation laws, and 
only the one case m = yields a nonlocal symmetry (previously found in 
[15115] ). 



Nevertheless, these potential systems may be very useful for finding new exact 
solutions [4] of the equations (NLS), (dNLS), (dNLS-H), (mKdV-l,2,H), which 
we will pursue elsewhere. 
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